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The end of the galactic cosmic ray spectrum 

Todor Stanev'' 

'^Bartol Research Institute, Department of Physics and Astronomy, 
University of Delaware, Newark, DE 19716, U.S.A. 

We discuss the region of transition between galactic and extragalactic cosmic rays. The exact shapes and 
compositions of these two components contains information about important parameters of powerful astrophysical 
sources and the conditions in extragalactic space as well as for the cosmological evolution of the sources of high 
energy cosmic rays. Several types of experimental data, including the exact shape of the ultrahigh energy cosmic 
rays, their chemical composition and their anisotropy, and the fluxes of cosmogenic neutrinos have to be included 
in the solution of this problem. 



1. INTRODUCTION 

In view of the expected high experimental 
statistics from the Auger Observatory ^ the in- 
terest in the transition from galactic to extra- 
galactic cosmic rays has increased significantly. 
It is now a standard belief that above some high 
energy of order 10^^ eV all observed cosmic rays 
should come from extragalactic sources because 
they cannot be contained in the Galaxy long 
enough to be accelerated For the rest of this 
paper we shall assume that the Ultra-High En- 
ergy Cosmic Rays (UHECR) are accelerated at 
powerful astrophysical objects. 

Extragalactic cosmic rays would lose energy 
in propagation from their sources to us if their 
sources are isotropically and homogeneously dis- 
tributed in the Universe. The main energy loss 
process is the photoproduction interaction in the 
microwave background (MBR) that causes the 
GZK effect, the steepening of the cosmic ray spec- 
trum above 6x 10^^ eV Several fits of the ex- 
isting experimental data have been published in 
recent years that derived different values of the 
most important astrophysical and cosmological 
parameters: the acceleration (injection in terms 
of cosmic ray propagation) energy spectrum of 
these particles and the maximum acceleration en- 
ergy, their chemical composition and the cosmo- 
logical evolution of their astrophysical sources. In 
the assumption that most of UHECR are pro- 



tons injection spectra as different as E"^ '' ^ and 
E~^-° and cosmological evolution of the type 
(1 + z)™ with m values from to 3 have been 
obtained. 

Other attempts |6I7I8| have assumed that ex- 
tragalactic cosmic rays have at their sources the 
same mixed chemical composition as the low en- 
ergy galactic cosmic rays. In such a case the 
main energy loss process is the disintegration of 
the heavy nuclei mostly in interactions in MBR. 
Hadronic interactions become important only af- 
ter the energy per nucleon exceeds the photopro- 
duction threshold. The fits of the observed cos- 
mic rays spectrum under this assumption gives 



an intermediate E 



injection spectrum. 



In all these attempts the fits of the observa- 
tions show the end of the galactic cosmic rays 
spectrum which is obtained by subtraction of the 
propagated extragalactic spectrum from the ex- 
perimentally observed one. This process gives 
some limits of the astrophysical parameters jH] 
when the subtraction gives unphysical negative 
values. 

An important and interesting question is what 
such models predict for the chemical composi- 
tion of UHECR, which for our purposes we shall 
define as cosmic rays of energy above 10^^ eV. 
The assumption that extragalactic cosmic rays 
are protons obviously lead to a composition that 
becomes very light at GZK and super-GZK en- 
ergies. Under the assumption that extragalactic 
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cosmic rays have a mixed composition the con- 
clusions are more comphcated and the question 
of the highest acceleration energy is more impor- 
tant. 

In terms of high energy astrophysics the loca- 
tion of the transition region can be used to es- 
timate the strength of the cosmic ray sources in 
different types of astrophysical objects, from spi- 
ral galaxies like our own to Active Galactic Nu- 
clei (AGN). In our current understanding the ac- 
celeration spectra depend on the strength of as- 
trophysical shocks and the maximum acceleration 
energy is a function of the magnetic field strength 
B and fluctuations SB/B. Other important pa- 
rameters are the cosmological evolution of the 
UHECR sources and the strength of the extra- 
galactic magnetic fields around the sources and 
the average magnetic field in the cosmologically 
nearby Universe. 

We shall discuss several types of data and their 
predictions in different models: UHECR spec- 
tra and chemical composition, the anisotropy in 
this energy range, and the production of cosmo- 
genic neutrinos by extragalactic cosmic rays. 

2. UHECR ENERGY SPECTRA 

Figure ^ compares two different fits of the ex- 
tragalactic cosmic ray spectra in the assump- 
tion that they are purely protons and that 
their differential acceleration spectrum is a power 
law E~'^~^^. The experimental data are from 
AGASA ^ and HiRes ^ and are normalized 
to each other at 10^^ eV. Since we are now only 
interested in the shape of the spectrum the exact 
differential fiux at the normalization points is not 
important. It is obvious that the two experimen- 
tal measurements agree well with each other on 
the shape of the energy spectrum with exception 
of the AGASA events above 10^° eV. The most 
recent analysis of the AGASA data, presented at 
this conference decreases the energy assign- 
ment of the AGASA data by 10-15% and makes 
the spectrum closer to that of HiRes. 

Fit a ^ derives an injection spectrum with 7 = 
1.7. The dip at about 10^^ eV is due to the transi- 
tion of proton energy loss to Bethc-Hcitler e+e^ 
pair production to purely adiabatic loss as pre- 



dicted by Berezinsky&Grigorieva 17 . The model 
does not need any contribution from galactic cos- 
mic rays to describe the observed cosmic ray spec- 
trum down to 10^^ eV. There is also no need for 
cosmological evolution of the extragalactic cosmic 
ray sources, although some source evolution can 
be accommodated with a slight change of 7. The 
model predicts purely proton composition of the 
extragalactic cosmic rays and does not work as 
well as shown in Fig. ^ if more than about 10% 
of the cosmic rays at the source are nuclei heav- 
ier than H. Fit b [1] does need contribution from 
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Figure 1. Comparison of two fits of the UHECR 
spectrum. Model a is from Ref. 4 and model b 
is from Ref. 0. See text for details. 

galactic sources with E~^-^'~' spectrum that ex- 
tends well above 10^^ eV. The extragalactic con- 
tribution is shown for two different cosmological 
evolutions with to = 3 (as used in Ref. 5 ) and 
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4 (upper edge of the shaded area). The influ- 
ence of the cosmological evolution on the cosmic 
ray propagation is modest because no protons 
injected at redshifts z higher than 0.4 arrive at 
Earth with energy above 10^^ eV independently 
of their initial energy. 

Obviously these two models predict very dif- 
ferently the end of the galactic cosmic ray spec- 
trum. In model a the galactic cosmic ray sources 
do not need to accelerate particles above 10^^ eV. 
In model h they should be able to reach energies 
higher by one and a half orders of magnitude. 
This would affect very strongly the expected cos- 
mic ray composition, as it will be discussed in the 
next Section. 

The assumption that extragalactic cosmic rays 
have a mixed composition at acceleration gives 
somewhat intermediate results for the injection 
spectrum of UHECR j6l7l8j . The spectrum that 
fits the observation best has a values between 1.2 
and 1.4. The chemical composition of cosmic rays 
at Earth is also different and obviously depends 
on the source composition. 

3. CHEMICAL COMPOSITION 

If the assumption that extragalactic cosmic 
rays are of mixed composition similar to that of 
low energy galactic cosmic rays is correct, that 
would change our perception for the energy de- 
pendence of the cosmic ray chemical composition. 
Before examining it let me first state that we are 
discussing the composition in terms of total en- 
ergy per nucleus, not the composition in terms of 
energy per nucleon as it is usually presented at 
energies between 1 and 100 GeV/nucleon. 

Since we believe that galactic cosmic rays are 
accelerated in stochastic processes at astrophysi- 
cal shocks (as they are in all models except in the 
Dar&DeRujula cannonball model |15ll()j ) we ex- 
pect cosmic ray nuclei to have maximum accelera- 
tion energy proportional to their rigidity, i.e. mo- 
mentum/charge. Similar dependence comes out if 
the knee of the cosmic ray spectrum is due to leak- 
age out of the Galaxy. The expectation is then 
that as soon as the maximum acceleration energy 
for hydrogen is approached the chemical 

composition starts becoming heavier - the max- 



imum acceleration energy for He E^%^^ is higher 
by a factor of 2 (factor of 4 in the Dar&DeRujula 
model). 

At the approach of E^^^ galactic cosmic rays 
become purely iron nuclei and this is supposed 
to be the cosmic ray composition when the tran- 
sition from galactic to extragalactic sources be- 
gins to happen. The measurements of the Kas- 
cade experiment |17) heavily support such behav- 
ior. Figure |21 shows the predicted composition 
as a function of the total energy per nucleus in 
the transition region in the classical cosmic ray 
units of < In A > . This is very appropriate when 
the detection is through air showers with logarith- 
mic sensitivity to both energy and composition. 
Model a presents the easiest case to explain. All 
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Figure 2. Cosmic ray chemical composition as a 
function of the total energy per nucleus in the 
three models discussed in the previous section. 
The mixed composition model data are taken 
Ref. 13 . The error bars assumed are logarithmic 
and probably lower than the true experimental 
ones. 



composition changes happen below 10 eV. We 
do not plot the composition below that energy 
as it is not very well defined. The results for the 
other two models are presented using our own un- 
derstanding of them. For model b we assume that 
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at 10^^'^ eV all cosmic rays are iron and assign 
< InA > value of 4. At higher energy we use 
the fraction of galactic cosmic rays to iron nuclei 
and thus calculate the corresponding < InA > 
value. For the mixed composition we use Fig. 3 
of Ref. 7 where the cosmic ray composition at 
Earth is plotted for 7=1.3 and E^ax given as 
1020-5 eV. 

Surprisingly the difference in the energy depen- 
dence of < InA > of model b and of the mixed 
composition model is not that big, although in 
the model b case the composition only consists 
of a combination of Fe and H, and in the mixed 
composition case we have five groups of nuclei. 
Distinguishing between these two models depends 
on the experimental sensitivity. If the composi- 
tion measurement uses the shower depth of max- 
imum Xmax as a function of the total energy of 
the shower, as done in fluorescence experiments, 
the error bars shown in Fig. [3 are approximately 
correct for SXmax of about 50 g/cnP. A better 
sensitivity should be able to distinguish between 
a two component composition of model b and the 
five component mixed composition model. 

Model a gives a very different picture, at least 
above 10^* eV, where the composition is purely 
Hydrogen. It is definitely distinguishable from 
the other two models. The prediction that the 
composition does not change above 10^* eV and 
is very light is supported by the HiRes X^ax 
measurement ^Hl- What is the exact meaning of 
light composition is not known because of the dif- 
ferences between the hadronic interaction models 
used for data analysis. 

On the other hand, other experiments support 
a much milder energy dependence of the cosmic 
ray chemical composition, that is more in line 
with the prediction of of models b and that of 
mixed extragalactic cosmic ray composition. 

4. ANISOTROPY 

Low energy cosmic rays diffuse in the magnetic 
fields of the Galaxy and lose memory of the lo- 
cation of their sources. The anisotropics are very 
small, well below 1% and are difficult to measure. 
The measurement of a small anisotropy with air 
shower experiments requires an exact knowledge 



of the experimental acceptance and of the lifetime 
of the shower array. 

At higher energy things are supposed to change 
as the cosmic ray rigidity increases and their diffu- 
sion in the Galaxy becomes faster. The first ques- 
tion we still cannot answer is of the rigidity (en- 
ergy) dependence of the cosmic ray diffusion co- 
efficient. From the ratio of secondary to primary 
cosmic rays an energy dependence of £;-0-5-0-6 
is derived. If this dependence is extended by 
seven orders of magnitude then galactic UHE pro- 
tons should show very strong anisotropy that is 
not observed. The two spots with increased cos- 
mic ray flux identified by the AGASA group ^H] 
(one in the general direction of the Galactic cen- 
ter and one in the direction of Orion) are not 
confirmed by other experiments. Theoretically 
we expect Kolmogorov turbulence in the Galaxy, 
that should under some circumstances give E^^^ 
energy dependence. 

One way to explain the low anisotropy in the 
beginning of the transition region is to assume 
that all galactic cosmic rays are heavy nuclei. If 
they were indeed heavy nuclei the average parti- 
cle rigidity would significantly decrease and would 
maintain the high isotropy of the galactic cos- 
mic rays. In the transition region, however, the 
cosmic ray chemical composition becomes lighter 
and correspondingly we expect to see some degree 
of anisotropy. Figure 13 shows the location of the 
sources of cosmic ray protons that reach the Earth 
in a propagation calculation. The source distri- 
bution is assumed to be Rqq where Roc is the 
galactocentric distance. There are no sources in- 
side the 4 kpc circle. Cosmic rays are isotropically 
injected in the galactic plane on a E~'^-'^ spectrum 
with energies between 4x10"'^^ and 2x10^^ eV. 
The magnetic field model is BSS with random 
field and a magnetic dipole in the center of the 
Galaxy. The Earth is approximated by a sphere 
of radius 1 kpc to speed the calculation. 

The figure shows the general behavior of pro- 
tons in this energy range. They do not diffuse, 
rather follow the magnetic field lines. Because of 
the magnetic turbulence the protons occasionally 
jump from one magnetic field line to another and 
may change the direction of their propagation de- 
pending on the pitch angle. 
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Figure 3. Location of sources from which isotrop- 
ically emitted protons in the galactic plane reach 
the Earth. See text for the details of the calcula- 
tion. 

When these protons arrive at Earth their ar- 
rival direction does not point at their source, they 
have rather a distribution that peaks in the direc- 
tion of the local field line. Since we expect higher 
source density inside the solar circle, the expected 
anisotropy may peak in the general direction of 
Orion. Exercises like that one do create relatively 
strong anisotropy, much higher than observed. 

One may complain that in Fig.|3|we are looking 
at particles that are accelerated in the galactic 
plane while in the case of extragalactic protons 
we will be dealing with particles that isotropi- 
cally enter into the Galaxy and may not create 
any anisotropy. The truth is that calculations like 
the one that generated Fig. O show that protons 
of energy above 10^^ eV emitted in the galactic 
plane leave the Galaxy much easier in the direc- 
tion of the galactic poles that after propagation 
parallel to the galactic plane. In the entrance 
of the Galaxy we would expect the same effect 
- more extragalactic protons should arrive from 
absolute high galactic latitude |6| than from low 
one. 

This would not be the classical anisotropy that 
surveys the arrival direction as a function of the 



galactic longitude, but it should still be significant 
effect. The details of the effect depend strongly 
on the galactic magnetic field model [201 and on 
details of the calculation such as the dimension 
of the 'Earth' in the calculation and of the prop- 
agation step size. 

I am convinced that in the case of high experi- 
mental statistics the anisotropy in the transition 
region between galactic and extragalactic cosmic 
rays deserves a careful study. A part of it is the- 
oretical. We should propagate particles of dif- 
ferent rigidity in the Galaxy and attempt to un- 
derstand their general behavior. Analytic solu- 
tions of the diffusion equation are not any more 
suitable for this problem. Particle trajectory has 
to be numerically solved, most likely in a Monte 
Carlo fashion and in detailed enough magnetic 
field models. Such models should be tested to 
match analytic calculations when applied to the 
appropriate simplified models. 

5. COSMOGENIC NEUTRINOS 

Cosmogenic neutrinos are produced in the same 
photoproduction interactions of the UHECR pro- 
tons that create the GZK effect. They were first 
proposed by Berezinsky&Zatsepin in 1969 [TU| 
and were the subject of many calculations after- 
ward. Cosmogenic neutrinos are often considered 
a 'guaranteed source' of ultrahigh energy neutri- 
nos. They are indeed guaranteed, since we know 
UHECR exist, but their flux is unknown. 

An essential quality of neutrinos is that they 
have a low interaction cross section. This makes 
neutrino detection a difficult problem that re- 
quires huge detectors of at least km^ scale. Such 
detectors are now in the stage of planning [23 and 
construction 22 23 . The question of the cosmo- 
genic neutrino flux and of its relation to other 
data from ultra-high astrophysics experiment is 
very timely. 

The low interaction cross section is not only 
a deficiency. While protons emitted at redshifts 
z exceeding 0.4 do not reach us with energy 
above 10^^ eV the cosmogenic neutrino produc- 
tion peaks at redshifts exceeding 2 for (1 -t- z)^ 
cosmological evolution of the UHECR sources. 
This is the main link between the extragalactic 
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UHECR and cosmogenic neutrinos. The detec- 
tion of cosmogenic neutrinos may help the degen- 
eracy in modehng of the extragalactic cosmic rays 
spectra shown in Fig. ^ 

The point is that models that use flat cosmic 
ray injection spectrum and require strong cosmo- 
logical evolution of the UHECR sources, such as 
model b would produce significantly more cos- 
mogenic neutrinos than steep injection spectrum 
models with no cosmological evolution. Such a 
comparison between the two models is shown in 
Figure 01 The figure shows cosmogenic neutrino 
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Figure 4. Fluxes of cosmogenic neutrinos gen- 
erated in proton interactions in the microwave 
background in the two models of the UHECR 



spectrum. Squares show the sum of 
the circles - the sum of 4- . 



and 



fluxes generated by proton interactions in MBR 
with the injection spectrum and cosmological evo- 
lution of the two models as indicated in the two 



panels. The injection spectra of the two models 
are normalized to the UHECR differential flux at 
10^^ eV. Before discussing the magnitude of the 
two fluxes I want to attract your attention to the 
flux of electron neutrinos and antineutrinos that 
exhibits two maxima separated by about two or- 
ders of magnitude of energy. The higher energy 
one is due to the t: i^^, + D^, + decay final 
state and consists mostly of electron neutrinos. 
It peaks at the same energy where /i^ and do. 
The lower energy maximum is due to electron an- 
tineutrinos from neutron decay. Comparison be- 
tween the neutron interaction length in MBR and 
their decay length show that all neutrons of en- 
ergy below lO^"'^ eV decay rather than interact. 

Model b generates much higher cosmogenic 
neutrino fiuxes than model a because of two rea- 
sons that contribute roughly the same increase of 
the cosmological neutrinos. Firstly, it uses much 
fiatter injection spectrum E^'^'^ which means 
equal amount of energy per decade. It thus con- 
tains many more particles above the photopro- 
duction interaction threshold with in the MBR 
is about 3x1019 eV and, of course, decreases as 
{l + z)-\ 

The other reason is that model b employs a 
strong cosmological evolution of the cosmic ray 
sources. This increases by the number of par- 
ticles injected at redshift of 2, but in addition 
increases by a factor of three the number of par- 
ticles above the interaction threshold. This way 
the total neutron flux is increased by the cosmo- 
logical evolution of the sources by a factor of flve. 

The difference in the peak values of the cosmo- 
logical neutrinos generated by the two models is 
more than one order of magnitude. In practical 
terms this means that model b generates fluxes 
that are in principle detectable by IceCube at the 
rate of roughly less than one event per year, while 
model a generates undetectable fluxes of cosmo- 
genic neutrinos in km^ detectors. 

Ice and water neutrino detectors are generally 
not very suitable for cosmogenic neutrino detec- 
tion. Much better strategies for these UHE neu- 
trinos are the radio and acoustic detectors that 
have very high detection threshold, but also will 
have higher effective volume. The other option 
are giant air shower arrays such as Auger, that 
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can reach effective volume of 30 km^ and, with 
sufficiently low threshold (10^^ eV) could detect 
several events per year. 

Cosmogenic neutrinos are also generated in the 
mixed composition scenario (25i26) . Since the 
major energy loss process is the disintegration of 
heavy nuclei, the main flux component of Ve from 
neutron decay. The Ve. flux exceeds by a factor 
of 5 the sum of the neutrino fluxes of all other 
flavors in Ref. The absolute magnitude de- 

pends again mainly on the injection spectrum and 
the cosmological evolution of the sources. 
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5.1. Production in the infrared/optical 
background 

The difference between different models is 
somewhat decreased when photon flelds differ- 
ent from MBR are considered. Several calcu- 
lations that include the infrared/optical back- 
ground (IRB) have been performed. The number 
density of IRB is of order 1 and varying by about 
a factor of 2 in different models, but its energy 
spectrum extends to energies as high as and ex- 
ceeding 1 eV. What that means is that extragalac- 
tic particles of much lower energy would interact 
in IRB and will generate neutrinos. Reference [77) 
shows the neutrino yields generated by protons of 
energy as low as 10^^ eV. This yield is small but 
has to be weighted by the much higher number 
of particles at that energy - almost 1,000 higher 
that that of 3x10^^ eV for a flat E~'^ injection 
spectrum. A more recent calculation that also 
includes UV photons 28 employs interactions of 
10"'^^ eV nucleons that have to be weighted still 
higher. Figure shows the v^^ ■\- cosmogenic 
neutrino fluxes from interactions in the MBR and 
IRB 1231 by the two models. The difference in 
the peak values is now slightly less than a factor 
of three. The steep injection spectrum model a 
has significantly more protons that interact in the 
IRB and for that reason partially compensates for 
the lack of cosmological evolution of the sources. 
On the other hand, since the average neutrino 
energy is substantially lower and the neutrino in- 
teraction cross section increases with energy, the 
event rates produced by the two models would 
not change much. If one uses steep injection spec- 
trum combined with an appropriate cosmological 



Figure 5. Fluxes of cosmogenic neutrinos gen- 
erated in the MBR and IRB. Only the sum of 
+ are shown. 

evolution of the sources that inclusion of the IRB 
target makes a big difference. 

6. DISCUSSION 

The transition region between galactic and ex- 
tragalactic cosmic rays is currently not very well 
defined and studied. It has a significant astro- 
physical importance, because the energy content 
of the extragalactic cosmic rays, as well as their 
injection spectrum and composition, will define 
much better the type of their sources. In addi- 
tion, a derivation of the cosmological evolution 
of the sources will not only restrict the number of 
available scenarios, but also provide an additional 
measure of the cosmological evolution of powerful 
astrophysical objects. 

This will not be easy to do even with the ex- 
pected much higher statistics and detection preci- 
sion of the Auger Observatory. We should employ 
all available types of measurements to solve this 
problem as well as it is necessary. These mea- 
surements include the cosmic ray composition in 
the transition region and above it, the cosmic ray 
anisotropics, and the possible detection of cosmo- 
genic neutrinos. 

The first experimental result on the shape 
of this transition presented by the Fly's Eye 
group 1221 used the simultaneous change of the 
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cosmic ray composition and energy spectrum 
shape. We have to do the same now - to com- 
bine aU observational data, look for consistency 
between different sets, and with models of the ex- 
tragalactic cosmic rays origin. 

The possible detection of cosmogenic neutrinos 
would be a powerful test if we succeed in collect- 
ing a reasonable statistics of such events. It is un- 
likely this will happen very soon, but the expand- 
ing efforts for designing and building detectors for 
ultrahigh energy neutrinos are very encouraging. 
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